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Abstract: Much attention has been given to dark matter explanations of the PAMELA 
positron fraction and Fermi electronic excesses. For those theories with a TeV-scale WIMP 
annihilating through a light force-carrier, the associated Sommerfeld enhancement provides 
a natural explanation of the large boost factor needed to explain the signals, and the light 
force-carrier naturally gives rise to hard cosmic ray spectra without excess 7r°-gamma rays 
or anti-protons. The Sommerfeld enhancement of the annihilation rate, which at low rela- 
tive velocities v re i scales as l/v re i, relies on the comparatively low velocity dispersion of the 
dark matter particles in the smooth halo. Dark matter substructures in which the velocity 
dispersion is smaller than in the smooth halo have even larger annihilation rates. N-body 
simulations containing only dark matter predict the existence of such structures, for example 
subhalos and caustics, and the effects of these substructures on dark matter indirect detection 
signals have been studied extensively. The addition of baryons into cosmological simulations 
of disk-dominated galaxies gives rise to an additional substructure component, a dark disk. 
The disk has a lower velocity dispersion than the spherical halo component by a factor ~ 6, 
so the contributions to dark matter signals from the disk can be more significant in Som- 
merfeld models than for WIMPs without such low-velocity ehancements. We consider the 
consequences of a dark disk on the observed signals of e + e~ , pp and 7-rays as measured by 
Fermi and PAMELA in models where the WIMP annihilations are into a light boson. We find 
that both the PAMELA and Fermi results are easily accomodated by scenarios in which a 
disk signal is included with the standard spherical halo signal. If contributions from the dark 
disk are important, limits from extrapolations to the center of the galaxy contain significant 
uncertainties beyond those from the spherical halo profile alone. 
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1. Introduction 



Much excitement has been generated over the past year or so by the release of results from 
several cosmic ray experiments. PAMELA [10, |26|, [| has shown a sharp upturn in the 
positron fraction above 10 GeV. Fermi ^ has published results indicating an excess in the 
total electronic spectrum e + + e~ in the 100 GeV to 1 TeV range, while HESS has shown 
that the spectrum falls off sharply above 1 TeV [[lj, Thermal WIMPs would naturally 
predict such features in the electron and positron measurements. However, conventional 
WIMPs annihilate far too little to yield these large signals, which require boost factors of 
0(100) above the thermal cross-section of 3 x 10~ 26 cm 3 s~ 1 [§l], |l|, ||, f|| 0. Moreover, 
the absence of an excess in antiprotons up to 100 GeV in the PAMELA || data essentially 
excludes models of dark matter which would produce the e + e~ excess through annihilations 
to standard model modes, since such models also have significant hadronic branching fractions 

An appealing possibility is for the dark matter (DM) to annihilate into a new, light 
(m < 1 GeV) boson p^]. This would result in a hard lepton spectrum and would kinematically 
prevent or suppress p production. Such annihilation modes easily fit the PAMELA and Fermi 
electron results [B9L B2|, while also producing an interesting synchrotron "microwave haze" 



signal, which has been observed by WMAP [56, 57, 48], and additionally yielding a significant 
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7-ray signal |49j from inverse Compton and final state radiation that may be observable at 
Fermi. A key consequence of including a new force in the dark sector is that at low velocities 
an annihilation rate much larger than that expected for a thermal WIMP can result from 
the Sommerfeld enhancement [[17], 76, 79] or from capture into bound "WIMPonium" states 



39, 90, pfl]. Sommerfeld-enhanced annihilation rates increase like l/v re i at small v re i, 



until the non-zero mass of the mediator eft cuts off the enhancement. As argued in 17], an 
important feature of such low- velocity enhancements is the added significance of substructure. 
Because the velocity dispersion in the dark disk can be significantly lower than in the smooth 
halo, up to ~ 6 times smaller, the annihilation cross-section and, therefore, annihilation rate, 
can be up to ~ 6 larger when the enhancement scales as \ jv re \. Add to this effect the increase 
in the annihilation rate due to the fact that the DM density in substructure is larger than in 
the smooth halo, and the role of substructure can indeed be substantial. 

The spectra of positrons and electrons can be greatly affected by the presence of sub- 
structure. Because distant [d ^ 1 kpc) positrons and electrons tend to lose their energy before 
reaching Earth, the high energy electronic signals, such as those measured by PAMELA and 
Fermi, are dominated by nearby sources. If there is nearby substructure, the observed spectra 
can be modified |3(|. Quantifying the effects, even for standard WIMPs, is a great challenge; 
the Via Lactea N-body simulation [^, [ft], ^] argues for a 1% chance of a subhalo giving a 
boost of 10 for conventional WIMPs, while the Aquarius simulation sees lower boosts [10C]. 
For WIMP scenarios with a Sommerfeld enhancement, quantifying the overall effect is more 
difficult. The Sommerfeld enhancement saturates at some velocity, below which the 1/v be- 
havior is no longer valid. If the saturation occurs at v ~ 200 km/s, as suggested in f99|| , 
then the annihilation cross-section is the same in the smooth DM halo, which has a velocity 
dispersion of ~ 220 km/s, as in substructure with lower velocity dispersion. If the saturation 
occurs at much lower velocities, the annihilation cross-section in substructure with the lowest 
velocity dispersion can be substantially larger than in the smooth halo. Since it is the smallest 
halos that have the lowest velocity dispersion, the smallest substructure can play a dominant 
role. This results in significant sensitivity to difficult-to-resolve small scales, and addition- 
ally to the unknown mass of the mediator <j>, which cuts off the enhancement, and cutoffs in 
the primordial power spectrum. Other complications, such as resonances in the annihilation 
cross-section for specific parameters, can make quantifying the effects nearly impossible. 

Nonetheless, it is important to attempt to understand what effects substructure can have 
on the present signals, in the context of Sommerfeld-enhanced annihilation channels. In this 
paper, we consider the effects of the addition of a "dark disk", as proposed by |)T], |92|], to 
the smooth spherical halo (see also ]p6[ ). Cosmological simulations that model only the dark 
matter particles are neglecting the baryons that dominate the mass of the galaxy interior 
to our location at the Solar Circle at R@ = 8.5 kpc. The dominant mass component is the 
stellar disk, and it plays an important role in the accretion history of the galaxy. |9l], [92] ] 
have shown that dynamical friction between the galactic disk and satellite galaxies causes the 
satellites to be preferentially dragged into the plane of the disk. Tidal forces subsequently 
disrupt the satellites, and the baryons and dark matter are accreted into disk-like structures. 
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The resulting dark disk, which is kinematically similar to the stellar thick disk, has a velocity 
dispersion of ~ 30 — 90 km/s. 

The presence of a dark disk in our galaxy can have consequences for the observed cosmic 
ray (CR) spectra. 1 The disk does not have the relatively large increase in mass density at the 
Galactic Center that the spherical halo has, and it extends only ~ 1 — 2 kpc from the Galactic 
Plane, so it sources fewer distant dark matter annihilation products. As a result, the disk 
gives a spectral hardening (as viewed here on Earth) of electrons and positrons, which are 
subject to large energy losses during propagation. See Fig. |l]a in which we plot the ratio of 
the dark matter positron flux due to annihilations from a pure disk to that from annihilations 
in a pure spherical halo, Disk /^Einasto- The ratio increases with energy, indicating that the 
relative contribution to the highest energy positrons is larger for the disk, which has a larger 
relative dark matter density in the Galactic Plane than a spherical halo. 

The spectral hardening of the electrons and positrons is not accompanied by a corre- 
sponding hardening of the proton and antiprotons, as the latter are not subject to the same 
large energy losses during propagation. In this way, substructure can change the relative 
numbers of annihilation products of a given energy reaching our detectors on Earth. See 
Fig. |l|b. The ratio of high energy antiprotons to positrons for the Einasto profile is ~ 1/3 
that for a dark disk. This can be understood in the following way. Compared to positrons 
of the same energy, high energy anti-protons have energy losses that are smaller by a factor 
of (^-) 2 , while their diffusion lengths are equal in the ultra-relativistic limit. 2 Within an in- 
finitely large and homogeneous diffusion zone, antiprotons diffuse a distance ^ (on average) 
larger than that of positrons with the same injected and observed energies. Models of our 
Galaxy indicate a diffusion zone of half-width L between ~ 2 and ~ 10 kpc. Protons and 
antiprotons of E ~ 100 GeV measured by PAMELA propagate to us from the entire volume 
of the assumed diffusion zone. On the other hand, electrons and positrons with energies ~ 1 
TeV measured by Fermi originate from a region within ~ 1 — 2 kpc of Earth, assuming their 
injection energies are not significantly higher. Thus, if a significant component of DM anni- 
hilation comes from a dark disk profile with a width of ~ 1 — 2 kpc, the observed antiproton 
flux will be suppressed relative to the predicted flux from a spherical halo alone, while the 
observed positron flux, being significantly more "local" in origin, will not deviate greatly from 
that predicted in the spherical halo scenario. 

The flux of photons coming from dark matter is highly dependent on the spatial distri- 
bution of the dark matter in the Galaxy. As an example, consider the differences between the 
fluxes expected from dark matter annihilation as compared to dark matter decay. Since dark 
disk profiles have a much lower density at the center of the Galaxy than Einasto spherical 
profiles |53| (see Fig. |2|), a dark disk that gives an 0(1) contribution to the local e^ 1 fluxes 
at E > 500 GeV, contributes significantly less to 7-rays and synchrotron radiation in regions 
within the inner 2 kpc. Thus, using local fluxes to find DM annihilation rates results in 

1 For a discussion of the effects of a dark disk on neutrino signals, see 

2 Within the validity of approximating the Klein-Nishina cross-section gk-n with the Tompson cross-section 

0~T, &KN ~ Ot- 



- 3 - 




Figure 1: Left: The ratio of the local flux of positrons from DM annihilation for the Dark Disk profile 
to that for the spherical Einasto profile. Here we assume monoenergetic electrons from the annihilation 
XX ~> e + e~ , with m x = 100 GeV, m x = 300 GeV, and m x = 1.0 TeV. All ratios have been normalized 
to 1 at the highest energy. Right: The ratio of the local flux of positrons {dot-dashed) and antiprotons 
(dotted) from DM annihilation for the Dark Disk profile to that for the spherical Einasto profile for 
the annihilation channel \X ~> W + W~, with m x — 100 GeV, m x = 300 GeV, and m x = 1.0 TeV. 
Also, we show the ratio of ratios (§p/$ e +) Disk /(&p/& e +) Einasto (solid). 



smaller fluxes of 7-rays and synchrotron radiation from the center of the Galaxy for dark disks. 
Because of this, a significant dark disk component in models with a Sommerfeld enhancement 
alleviates some the constraints coming from gamma ray measurements of the Galactic Center 
and Galactic Ridge regions. As a consequence, while the Galactic Center places the tightest 
constraints on the annihilation rate for cusped, spherical profiles, the Galactic high latitudes 
(specifically, the isotropic diffuse flux there) usually places the tightest constraints on the 
annihilation rate for the disk profiles. 



2. Definitions and Parameters 



2.1 Dark Matter Halo Models 



In the ACDM cosmological simulations of [ 91 1 , the formation of the stellar thick disk through 
accretion can be accompanied by the formation of a dark matter disk. The fits of the dark 
disk (DD) component can be described by an exponential function and parametrized by half 
mass scale lengths. We consider a disk component based on the fits of [^IJ, and assume a 
cylindrically symmetric dark disk density profile of the form 



p(R, z) = p exp 



1.68 (R Q - R) 



R1/2 





0.693 \z\~ 


exp 






z l/2 . 



(2.1) 



- 4 - 




o.o L -I o.o L - 

5 10 15 20 1 2 3 4 5 



cylindrical radial distance R (kpc) at z-0 distance z (kpc) at R-8.5 kpc 

Figure 2: Left: Radial dependence of DM density for the Einasto profile and the Dark Disk profile 
at z = 0. Right: z dependence of DM density for the Einasto profile and the Dark Disk profile at 

T cylindrical "0 • 



where -R1/2 = 11-7 kpc and ziy 2 = 1-5 kpc are the half mass scale lengths in the Galactic 
Plane and perpendicular to the Galactic Plane, respectively, and Rq = 8.5 kpc. Here R is 
the cylindrical radial coordinate. 

For the spherically symmetric halo model (SH) we use the Einasto profile |53|] following 
Merritt et al. 



p(r) = po exp 



2 (r a -R% 



a 



-2 



(2.2) 



with a = 0.17 and r_2 = 25 kpc. r is the spherical radial distance from the Galactic Center 
(GC). We take the value of the local density for both the spherical and disk components to 
be po = 0.4 GeVcm~ 3 [27]. (See discussion below.) In Fig. || we show the dependence of 



these DM profiles on R for z = 0, i.e. in the Galactic Plane, and on z for R = 8.5 kpc, i.e. 
perpendicular to the Galactic Plane at our location. 

According to |)l]], values for the ratio of the local DM density in the dark disk to the 
local DM density in the spherical halo Po dd /po S h ran S e from 0.2 up to more than 1.5. The 
highest values of po DD / Po S h occur i n simulations when the thick stellar disk has a higher mass 
density than the thin stellar disk. This can result from an increase in the mass included in 
the stars of the thick disk, or from a thin disk heated by very massive, high-redshift mergers, 
which results in the thick disk getting populated by thin disk stars. Another, though less 
likely, cause is the occurrence of multiple prograte and low inclination mergers [91]. 



Constraints have been placed on the total local dark matter density from estimations of 
the local dynamical mass density based on measurements of proper motions and parallaxes 
of stars. Using the Hipparcos data, Q] derived the local dynamical mass density to be 
0.102 ± 0.010 M Q pc~ 3 (95% C.L.). Their estimate for the local density of the visible matter 
is 0.094 ± 0.017 M Q pc -3 (67% C.L.). Assuming the mean values of these densities, we 
get 0.008 Mqpc~ 3 , or equivalently 0.3 GeVcm~ 3 , for the local dark matter density with 
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an error that is larger than the mean value, thus allowing for a dark disk with pq l 



JDD 

-3 



Po SH , where Po SH = 0.4 GeVcm . [44| derived a smaller local dynamical mass density of 



0.076 ± 0.015 M Q pc 3 using the Hipparcos data, but values of 0.3 — 0.4 GeVcm 3 are allowed 



within their uncertainties. Later analysis of the Hipparcos data by [64] and [72] agree with 
the constraints imposed by |j3| . Recently |27]] suggested a value of 0.4 GeVcm -3 as the best 
estimate of the local DM density. Also |95| calculated the value of the local DM density to 



be 0.43 ± 0.11 ± 0.10 GeVcm -3 , while |103j ] set an upper limit of 1 GeVcm -3 for the local 
DM density. 

2.2 Cosmic Ray Propagation 



We use GALPROP, developed by Moskalenko and Strong [101], for the propagation of all 
cosmic ray particle species. For our calculation of the Inverse Compton Scattering of elec- 
trons and positron we use the Interstellar Radiation Field model of [p9| . For consistency in 
our comparisons of the spectra from a dark disk and from a spherical halo, we assume the 
same energy losses, diffusion coefficient, and diffusion zone 3 . Our choices for the background 
models and propagation parameters give local CR spectra that agree with local cosmic ray 
measurements, including the proton, He, C, Fe, B/C, Be 10 /Be 9 , and sub-Fe/Fe spectra. Addi- 
tionally, for the annihilation channels with no DM contribution to antiprotons our background 
antiproton-over-proton ratio p/p agrees with the PAMELA p/p data. See Fig. [|a. 

In placing constraints on p fluxes from DM, we use a background antiproton ratio that 
is 15% smaller. The smaller p/p background ratio shown in Fig. |9]b is consistent with the 
current constraints from all CR data. For example, a reduction in the background p/p could 
be based on uncertainties in the interstellar medium (ISM) distribution of Hydrogen and 
nuclei, since antiprotons are produced by collisions of CR protons and nuclei with the ISM 
protons and nuclei. 

2.3 Annihilation Channels 

We consider two classes of annihilation scenarios, the "eXciting Dark Matter" (XDM) scenario 
[|l7], [}(], 58] in which annihilation proceeds through a single mediator (j), and annihilation 



proceeding through two mediators <3? and 4>. Additionally we consider annihilation directly 
into muons, XX A t+ / U_ ; an d annihilation into VF+I'F - , which then decay. We discuss the 
implications of the W + W~ channel for the PAMELA antiproton data. 
In the XDM scenario, we present results for annihilation 

i) into e + e~ through a scalar (or vector) (j) with mass 2m e < < 2m^ and a branching 
ratio (BR) of 0.9 4 (or 1), 

ii) into p + p~ through a scalar (or vector) <p with mass 2m ^ < < 2m n and BR=1, 
hi) into 7r + 7r - through a vector (f> with mass = 750 MeV (the form factor for <p — > 

vr+ + vr - peaks at ~ 750 MeV |52|), 



3 Refer to [g2| or |l02{ for the details of the diffusion parameters. 
4 We expect a ~ 10% decay of scalar <f> — ► 27 due to 1-loop diagrams. 
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iv) and into e + e~, p + p~ and ir + ir~ with relative BR's 1:1:2, through a vector <f) with 
mass ~ 650 MeV. 

For DM annihilation with subsequent cascade decay through two mediators, we show the 
case XX ~^ with subsequent decays $ — > 4>4> and 4> — > where 2m^ < < 2m % 

and m$ ~ 10m^, and additionally the similar cascade decay with <f> — >• e + e~. The spectra of 
the final e + e~ injected into the ISM from these cascade decays are given in Appendix A. 

2.4 Annihilation Rate and Boost Factor 

In the absence of the Sommerfeld enhancement, the total dark matter annihilation rate T ann 
from the spherical halo and dark disk is given by (for Major ana fermions) 

r ^(psh + Pdd\ 2 , | ,v , nQ , 

Tann = ~ \ a ann \V ). (2.3) 

2 \ m v / 



{&ann \ V \ ) mixed • 



If the Sommerfeld enhancement is present, Eq. 2.3 is modified to 

r 1 / PSH \ 2 / | I \ . 1 ( PDD\ 2 , | |> . ( PSH ■ PDD\ 

(2.4) 

(fannl^Dsi/ and (<7 ann |?;|)i5D are in general different, because the velocity dispersion of dark 
matter is different in the spherical and disk components of the halo. 5 

We define the boost factor (BF) as the ratio of the thermally averaged annihilation cross- 
section needed to fit a set of data, (cr ann \v\) fu, to 3 x 10~ 26 cm 3 s _1 , the expected thermally 
averaged annihilation cross-section for a WIMP with mass m x ~ 500 GeV, 

BF= ( ga ""H>/; (2.5) 
3 x 10" 26 cnrV 1 V ' 

3. Results 

3.1 Positrons and Electrons 



Sommerfeld Enhanced Annihilation Channels 

As mentioned earlier, the effect of a dark disk on the electronic signal is to produce a 
harder spectrum at high energies relative to that produced by the spherical Einasto profile, 
bee Fig. |a. This is a generic result for the dark disk profiles of and any spherical halo 
profile, and it arises because lower energy electrons and positrons, which propagate to Earth 
over larger distances on average, are probes of regions closer to the Galactic Center where 
the spherical halo profiles peak significantly more than the disk profiles. See Fig. ||a. This 
also explains why the boost factor needed to get equal fluxes of dark matter electrons and 
positrons at energies much less than m x is lower for the Einasto profile than for the dark 
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Figure 3: Positron fraction $ e + / ($ e + + $ e - ) (left) and total electronic flux weighted by 
E 3 (GeV 2 m~ 2 s _1 sr _1 ) (right), both as a function of energy. Results are shown for the spherical 
Einasto profile (dashed) and for a dark disk (solid). Backgrounds are shown in dot-dashed. Top: 
XDM fj+fj,- channel with m x = 2.5 TeV. Bottom: XDM tt+tt" channel with m x = 3.1 TeV. PAMELA 
postron fractio data are from || , while the total electronic flux measured by Fermi given in |Q, . 



disk, i.e. why the boost factors needed to fit PAMELA are lower for the Einasto profile than 
for the dark disk. 

In Figs. |3||8] we show the positron fraction $ e +/($ e + + $ e -) and the .E 3 -weighted total 
electronic flux (the sum of the fluxes of electrons and positrons) for a number of annihilation 
channels. We show the fits of the spectra to the PAMELA, Fermi, and HESS data separately 
for the dark disk and spherical halo scenarios. The locally measured spectrum, which is 
composed of both a disk and a spherical halo component, is expected to lie between the two 

5 For the mixed term (oc psn ■ Pdd), we take (a\v\) mixe d = {(T ann \v\)sB- 
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Figure 4: Cosmic ray signals as in Fig. ||. XDM e + e , fi + fi and 7r + 7r channels with relative BR's 
of 1:1:2. 



spectra shown. 

For the relatively hard annihilation channels, such as \X shown in Fig. [?], the 

DM component of the spectrum arising from the dark disk is too peaked between 400 and 
800 GeV to give really good fits to the Fermi data, though the fits for the Einasto profile for 
these channels are generally good. Those annihilation channels with much softer electron and 
positron spectra, for example XDM to muons or to pions only or some combination of these 
and XDM through a 2-step decay (see Figs. H||5|), give very good fits to the Fermi data for 
both the spherical and disk profiles. Therefore, any relative combination of disk and spherical 
halo contributions to the total e^ flux is allowed by the Fermi data. (See also Section |3.3| 
for more discussion on this). Both DD and SH profiles give good agreement with the sharp 
fall-off in the spectrum seen in the HESS data, which suggests a power law of E 



-3.9 ± 0.1 : 



111 



the range 700 - 3500 GeV §\A\. 

We note that for the flux of primary electrons we use an injection power-law of ~ E~ 2 ' 5 , 
which agrees with conventional assumptions about the primary electron spectrum in the GeV 
energy range [77, ||, 71]. We introduce a break in the spectrum at ~ 2 TeV in order to 
have agreement with the HESS data. Since the primary electrons composing the bulk of 
the background flux at TeV energies are ISM electrons accelerated from supernovae shocks, 
a break in the spectrum at this energy scale is natural. However, the energy at which the 
break occurs depends on the most recent and closest supernovae events and is not known 
theoretically; rather, it comes from fitting local CR data. 

The PAMELA positron fraction data is well fit for the Einasto profile by all of the XDM 
annihilation channels shown. The implication for the positron fraction of the harder DM 
e + e~ spectra for the disk profile is a steeper rise of this fraction. For the dark matter masses 
that fit the Fermi data, m x > 1.0 TeV, this can result in slightly worse fits to the PAMELA 



-9- 



data. See, for example, Fig. |3[ However, for low DM masses that agree with PAMELA but 
not Fermi, the dark disk can improve the PAMELA fits. 

Since the velocity dispersion of the particles in the dark disk traces that of the stars in 
the stellar thick disk, the local values are much lower than for the spherical Einasto profile. 
Typical local values of the velocity dispersion for the disk are ~ 30 km/s and for the spherical 
halo are ~ 220 km/s. Therefore, we expect that if the annihilation cross-section scales as 
~ 1/v, then in the disk it can be up to ~ 10 higher than that of the particles in the spherical 
halo. This is assuming, of course, that the Sommerfeld enhancement isn't saturated for 
velocity dispersions greater than that in the disk and that we are not near a resonance. (In 



Section |3.3| , we discuss the possibility, as suggested by 9S ] , that the annihilation cross-section 
is close to saturation for a v ~ 200 km/s.) In Fig. || we show an example of the PAMELA and 
Fermi spectra for a combination of disk and Einasto profiles for a disk annihilation rate that 
is 5 times the annihilation rate in the spherical halo. 

For a specific DM mass and annihilation channel, the ratio of the boost factor needed to 
fit the PAMELA data over the boost factor needed to fit the Fermi data, BFpam /BFp erm i, 
is larger for the dark disk profile than for the Einasto profile. For example, for the XDM 
7r + 7r~ channel BFpam / BFpermi = 3.2 for the dark disk, while for the Einasto profile the 
ratio is 1.9. Due to the spectral hardening of e^ of DD profiles compared to SH profiles, 
the ratio of boost factors B Fpam / B Fp erm i is expected to be higher for DD profiles than 
for SH profiles. In other words, if the fluxes of e + e~ from the spherical halo and from the 
dark disk are equal at 100 GeV, then the flux due to the spherical component will be greater 
for E < 100 GeV, while the flux for the dark disk component will be greater for E > 100 
GeV. Thus, a spectral hardening would allow annihilation channels that under-produce at the 
highest Fermi energies in a spherical halo to have improved fits to the data for a disk halo. 
On the other hand, the "flatness" of the Fermi data (in £' 3 -weighted plots) can be used to 
constrain the disk contribution to the overall flux. Additionally, WMAP data can put upper 
limits on the saturated Sommerfeld enhancement [99], so the contribution of the dark disk 
for various channels can be constrained in this way as well. See Section 3.3. 



Non-Sommerfeld Enhanced Channels 

The annihilation channels that proceed through one or more light mediators (see Figs. |||- 
^) require large boost factors to fit the PAMELA and Fermi data, and these can be explained 
by the Sommerfeld enhancement. Annihilation channels directly to leptons also require large 
boost factors, BF ~ 500 for the spherical halo, but in these scenarios there is no long range 
force to give rise to the Sommerfeld enhancement needed to motivate such a large boost factor. 
A Breit-Wigner resonance, as suggested by could account for the large annihilation rate 
needed at z = to produce the signals we observe today, while not impacting the thermal 
relic density at freeze-out. Also, WIMPs produced non-thermally, such as those of |5l|, [23[ |, 
can naturally have large annihilation cross-sections that give good fits to the PAMELA and 
Fermi data. For this reason, we consider here such non-Sommerfeld enhanced scenarios in 
the context of a dark disk. 
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Figure 5: Cosmic ray signals as in Fig. |3|. Top: XDM e + e~ through two steps with m$ ~ lOm^ and 
rrirj, < 2m^. Bottom: XDM through two steps with 2m^ < < 2m w . Two-step annihilations 

proceed through \X ~^ followed by the decays $ — > cjxf) and — > e + e~ or <\> — > 



In Fig. |?| we show the positron fraction and total electronic flux for dark matter annihi- 
lating monochromatically to muons, \X ~> - Good fits to both the PAMELA and Fermi 
data can be achieved with the spherical halo profile. As mentioned earlier, the disk profile 
gives an electronic spectrum that is too hard to fit well the Fermi data, though it does a good 
job reproducing the positron fraction. Additionally, the ratio BFpam / BFF erm i is 4.5 for the 
dark disk, which suggests that a dark disk dominant scenario does not accurately explain 
both signals simultaneously. 

While direct decay to taus, \X ~> t + t~~ , gives good fits to the Fermi and PAMELA data, 
the annihilation rate for the tau channel is highly constrained by the Fermi diffuse 7-ray data 
[ [3^ , |85|| and by the HESS measurements of gamma rays in the Galactic Ridge 



12, |2|. (The 



- 11 - 



(a) 



Halo + Disk, BF = 140 
Background 
PAMELA Data 




XDM e + e, 1 :1 Channel 
m v = 1.3 TeV 



E 

°> 
o 

o 100 



(b) 

HESS Systematic Error Band 

HESS Data 
— i — i Fermi High E Data 
--i---' Fermi Low E Data 

Halo + Disk, BF = 45 

Background 




XDM e e", [l fi" 1 :1 Channel 
m y = 1.3 TeV 



Energy (GeV) 



100 1000 
Energy (GeV) 



Figure 6: Cosmic ray signals as in Fig. g. XDM e + e~ and channels with eqnal BR's. Here a 

combination of dark disk and Einasto profile is assumed. The annihilation rate in the disk is taken to 
be 5 times the rate in the spherical halo. 
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Figure 7: Cosmic ray signals as in Fig. ||for annihilation of DM directly to for rn x = 1.5 TeV. 



t + t channel produces a high energy prompt component of gamma rays through the 7r°s 
created in the decay chain.) Annihilation directly into e + e~ gives a propagated spectrum 



e + e with spectral index of -2 and a hard cut-off at E = m x [29, 78|; such a spectrum does 
not fit the Fermi and HESS e~ + e + data well. 

3.2 Antiprotons 

Annihilation through W + W~ gives very poor fits to the Fermi and PAMELA data. Most 
of the electrons and positrons are produced in the decay chains of the hadrons and taus, so 
they have relatively soft spectra. The resulting abundance of low energy e + e~ pairs gives a 
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Figure 8: Cosmic ray signals as in Fig. | for DM annihilation to W + W~ with m x = 1.0 TeV. 
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Figure 9: Antiproton over proton ratio p/p as a function of energy. Left: Background for the 
propagation and diffusion parameters used in our calculations. Right: For DM annihilating through 
W + W~ with m x = 1.0 TeV for the Dark Disk profile (solid) and for the Einasto profile (dashed). Here 
the background (dot-dashed) is smaller than the standard background by 10%. Data of the antiproton 
to proton flux ratio are from B3]. 



positron fraction that is much too flat and an e + + e~ spectrum that doesn't simultaneously 
agree with the Fermi data around 20 GeV and the data around 1 TeV. In Fig. || we show the 
results for m x = 1 TeV, the DM mass that best fits the Fermi data. The fit to the Fermi 
data is significantly better for the dark disk profile than for the spherical halo. Additionally, 
the fit of the positron fraction to the PAMELA data has a modest improvement. Thus, 
for the W + W~ channel the favored scenario is one in which the dark disk is the dominant 
contribution to the local DM density. 
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An additional benefit of the dark disk for the W channel is that the constraint placed 
on the annihilation rate by the PAMELA antiproton data is somewhat relaxed. In Fig. [^b 
we show that a very similar antiproton ratio is obtained for the dark disk with BF = 150 
as for the Einasto profile with BF = 25. The factor of ~ 6 "improvement" is a reflection 
of the fact that the p's propagate to us with little energy loss from the entire diffusion zone, 
and the disk has a smaller average dark matter density throughout the diffusion zone, as 
discussed in Section [j]. This, of course, depends greatly on the relative widths of the diffusion 
zone and the dark disk. For diffusion zones with widths much larger than the scale width 
of the disk, the improvement becomes significant, i.e. factors are much larger than 1. For 
narrow diffusion zones such as the ~ 1 kpc suggested by |57]], the thickness of the disk is 
irrelevant, as it is thicker than the diffusion zone. If the dark disk is thicker than the diffusion 
zone, improvement can still be had if the average value of the DM density over the diffusion 
zone is less for the dark disk than for a spherical halo profile. We note that for m x = 1.0 
TeV annihilating through W + W~ in the dark disk only scenario, the BF needed to fit the 
PAMELA positron fraction is still ~ 4 times larger than that allowed by the PAMELA p/p 
ratio. For comparison, in the spherical halo only scenario the BF for the PAMELA positron 
fraction is ~ 10 times larger than that allowed by p/p. However, since the fits to the PAMELA 
data are so poor, the -BF's should not be taken too seriously. In summary, the discrepancy 
between the larger boost factors needed to fit the positron fraction data and the smaller 
boost factors needed to fit the antiproton fraction data can be reduced by including a disk 
component. Thus, a higher branching ratio to hadronic channels is allowed for a dominantly 
dark disk profile. 

3.3 Constraints on Boost Factors 



Constraints from the Fermi e + + e flux 

The combined flux from a spherical halo and a dark disk can be written as 1 



dN tota i dNprim dN sec BF S h (dN S H , D ,, 2 dN DD „,„dN comb \ 
-^ = ^E- + ^E- + JTTW\~^ -dE- + 2M ^E-)- (3 ' 1} 

Here M is the ratio of the local DM densities, Po DD / Po SH , R is the ratio of the thermally- 
averaged annihilation cross-sections, R = {o-\v\)dd/{&\v\)sh, and BFsh is the Sommerfeld 
enhancement in the spherical halo (i.e. for the velocity dispersion in the spherical halo). 

Using the e + + e~ flux as measured by Fermi, we can constrain M, R, and BFsh- Since 
the Fermi data constrain the total flux of Eq. 3.1, constraints on any of these parameters can 
be used as further input to constrain the rest of the parameters. We consider the following 
values for the parameters: 

< M < 1, POdd < p 0sH , p 0sH + p 0DD = 0.4 GeVcm" 3 . (3.2) 



The — jfp^ comes from the term oc pDD * Psh in Eq. 2.4, where we assume that the relative velocity 
dispersion between the particles in the DD and the SH is ~ 200 km/s. 
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We take R m ax = 100 to be the maximum ratio of the annihilation cross-section of DM 
particles in the disk halo to the spherical halo. This value could occur if we live close to a 
resonance, where the enhancement scales as l/v 2 , since the velocity dispersion in the spherical 
halo is ~ 220 km/s and in the dark disk halo is ~ 30 km/s. 7 This assumes that the Sommerfeld 
enhancement does not saturate at velocities above ~ 30 km/s. If saturation occurs for 
annihilation in the dark disk but not in the spherical halo, i.e. at a velocity between ~ 220 
and ~ 30 km/s, then R is the ratio of the saturated cross-section in the disk to the cross- 
section in the spherical halo. 

In our calculations we allow the normalization of the primary e~ flux to vary in the range 
1.48 - 2.04 x 10" 2 GeV _1 m -2 s _1 sr _1 at 20 GeV, and the power law index (dN/dE ~ E~ a ) 
to vary in the range 3.2 < a < 3.35. Since the secondary e s contribute negligibly to the 
total e + + e~ flux for E > 20 GeV, we hold them constant assuming dN sec /dE ~ E~ and 
a normalization of 6.7 x 10 -4 GeV m~ 2 s sr _1 at 20 GeV. 

In column 4 of Table [TJ we show the 95% C.L. values for R, the ratio of the annihilation 
cross-section of DM particles in the disk to the spherical halo, assuming that the local DM 
density receives equal contributions from the disk and the spherical halo. If the Fermi error- 
bars are correlated, the x 2 distribution test is not an exact measure to calculate confidence 
levels. Since we don't have any information apart from the errorbars, the mean values of the 
fluxes, and the energy binning, we consider each point as an independent normally-distributed 
variable with mean and standard deviation given by the Fermi data. For all cases except the 
XDM electron channel, the channel with the hardest DM e^ spectra, the annihilation rate in 
the disk can be at least 100 times the annihilation rate in the spherical halo, indicating that 
the Fermi data is consistent with the presence of a dark disk component in a non-saturated 
resonance scenario. In column 5 of Table p] we show the 95% C.L. for the maximum value 
of M, the ratio of the local DM densities, assuming that the annihilation rate in the disk 
is 100 times the annihilation rate in the spherical halo. The Fermi e + + e~ data does not 
constrain the ratio of the local DM densities within our allowed range of values (0 to 1), 
except for the XDM electron channel. The results of columns 4 and 5 in Table |l] suggest that 
for the XDM channels giving the best fits to Fermi and PAMELA, the dark disk may be a 
dominant contributor to the local flux of e^ coming from dark matter annihilation. In the 
latter case, the Sommerfeld enhancement in the spherical halo is typically smaller, by an order 
of magnitude, than the enhancement when the spherical halo is the dominant contribution. 

In the limiting case, Po DD /po SH "C 1, the dark disk contribution can be significant to the 
flux of e^ locally only if the masses of the dark matter m x and the mediator and the value 
of the coupling between the two A are such that for a velocity dispersion of ~ 30 km/s the 
annihilation cross-section has a resonance fT7| , |88| . The results shown in Table [l| indicate that 
XDM channels that produce e ± with a spectral index a > 2 at 300 GeV 8 , where ^ oc E~ a , 

7 Since we want to place an upper bound on the Sommerfeld enhancement in the dark disk, we use a value 
for the velocity dispersion that is on the low side, but within the range suggested by simulations Ell. Such a 
value is motivated in dominant thin disk scenarios. 

8 XDM to pions only, XDM 2-step to muons, and XDM to e , fj, and tt^ 1 with relative ratios of 1:1:2 
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Channel 


M x (TeV) 


BF SH 


Rmax{M = 1) 


M max (R = 100) 


BF\YMAP5 


XDM to e ± 


1.2 


10 - 160 


90 


0.9 


210 


XDM to /i ± 


2.5 


30 - 1100 


100 


1 


1300 


XDM to vr ± 


3.1 


40 - 2000 


100 


1 


1900 


XDM 1:1:2 


1.6 


7.0 - 380 


100 


1 


560 


XDM 2-step to e ± 


1.5 


2.0 - 250 


100 


1 


260 


XDM 2-step to /x ± 


3.8 


10 - 2100 


100 


1 


1900 




1.5 


220 - 1550 


1* 


1 (R = 1) 


750 



Table 1: Table of allowed values of the boost factor for the Einasto profile for fits to the Fermi data, 
BFsh, the maximum allowed values of R = {<j\v\)dd/ {cf\v\)sh assuming M = po DD /po S H = 1 (95% 
C.L.), the maximum allowed values of M for R = 100, and the boost factors BFwmap5 excluded by 
WMAP5 data at 95%C.L. *There is no Sommerfeld enhancement for \X ~ * P + P~ ■ 

are allowed within 95% C.L. to have a resonance in their cross-section. We clarify that in the 
cases of resonant Sommerfeld enhancement at z = for a velocity dispersion of 30 km/s, it 
is the dark disk that gives the dominant contribution to the local e ± flux, while the spherical 
halo is constrained to be subdominant. See Table |l[ At v re i ~ 200 km/s for m x ~ TeV, the 
Sommerfeld enhancement is O(10). 

Since the PAMELA positron fraction has only 6 data points at energies E e ± > 10 GeV, 
it constrains the annihilation rates less than the Fermi e + + e~ flux. Additionally, the fact 
that the secondary e + 's are still at least half of the total e + flux for E e + < 50 GeV allows for 
greater uncertainty in the necessary DM annihilation BF's. 

Constraints from WMAP5 data 

The high energy electrons and positrons produced by dark matter annihilation in the early 
universe heat and ionize the photon-baryon plasma and thus can give rise to an increase in the 
ionization fraction of the plasma after recombination. Perturbations in the ionization history 
around the time of last scattering (z ~ 10 3 ) can affect the temperature and polarization power 
spectra of the CMB. [|99| (see also |105f| ) calculated the effect of annihilating DM on the plasma 
and placed constraints on the maximum annihilation cross-section using the WMAP5 data. 

The upper bounds on the cross-section depend linearly on the mass of the DM candidate 
and are inversely proportional to the efficiency / with which the energy from DM annihilation 
gets converted into the energy of the e ± products. 9 In column 6 of Table [l] we present the 
upper bounds on the boost factor (assuming {a ann \v\} = 3 x 10~ 26 cm 3 s _1 ) for several XDM- 
type annihilation channels. For the XDM e^, XDM fi^, XDM ir^ channels, we used the mean 
efficiency / provided in Table 1 of [p9|[ . For XDM 1:1:2 the efficiency is a linear combination 
of the efficiencies for XDM e ± , XDM /i ± and XDM vr ± , while for XDM 2-step to e ± and 

9 The efficiency / depends on the annihilation channel, the redshift, and also has a weak dependence on the 
mass of the DM candidate. 
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2-step to , the efficiency is approximately equal to the efficiency of XDM e ± and XDM p^ , 
respectively. 

The values of BFwmap5 presented in the table are upper limits on the allowed BF from 
Sommerfeld enhancement at z ~ 1000. Since the velocity dispersion of DM in both the 
spherical halo and dark disk are expected to be larger now than at z ~ 1000 (resulting in 
smaller boost factors now), the values in Table [I] can be considered the maximum allowed 
values of the boost factors today (95 % C.L.). Current values of the annihilation cross-section 
would be equal to those at z ~ 1000 if we are in the saturated regime. 

A comparison of the BF's for the fits shown in Figs. with the allowed ranges of BF$h 
and R of Table [l] makes it clear that the CMB induced limits are low enough to constrain the 
contribution of the dark disk for the various channels. We remind the reader that the BF's 
given in Figs. ^-^ assume a local DM density of po = 0.4 GeVcm -3 , and that, as discussed 
in Section ^ the value could be different by a factor of ~2. If po = 0.6GeVcm -3 , the BF's 
needed to fit the Fermi or PAMELA data would be a factor of ~ 1/2 of those shown. Thus, 
a factor of a few discrepancy between the fitted BF and the CMB limit is allowed. 

For the XDM p^ channel, the Fermi data excludes neither the spherical Einasto halo 
nor the dark disk. The boost factor for the fit of the dark disk to the PAMELA data is ~ 5 
times bigger than the allowed value. 10 (Similar results hold for the XDM e channel whose 
fits to the data are not shown.) However, if the contribution to the local e^ flux of the dark 
disk and the spherical halo are comparable at ~ 500 GeV, the fluxes at PAMELA energies 
come mainly from the spherical halo. Thus, even if the dark disk gives an 0(1) contribution 
to the local flux at E ~ 10 3 GeV, those channels are not yet excluded by the CMB limits, 
but are in some tension. 

The XDM tt ± and XDM 1:1:2 channels are within CMB constraints, assuming the factor 
of 2 discrepancy between the fitted BF and the CMB limit mentioned earlier. In the com- 
bination channels, less of the available energy from DM annihilation is dumped into at 
z ~ 1000, so there is a smaller effect on the CMB temperature and polarization power spectra. 
For the XDM 2-step to e ± and 2-step to p ± , the PAMELA fits require BF's within a factor 
of a few of the CMB limits, thus those channels aren't excluded. We note that regardless of 
the specifics of the annihilation channel (and the needed DM mass and BF), the fits to the 
PAMELA and Fermi data constrain the spectrum of DM e ± to be similar in all cases, thus 
we expect the CMB limits to be similarly constraining in all cases. 



Recently |55|, 54] suggested stronger constraints on the Sommerfeld enhancement based 



on relic density calculations. In |55j the authors suggest that for m x = 1 TeV and = 250 
MeV, the maximum allowed Sommerfeld enhancement is 90. In general, they find that for 
m x > 1 TeV and ~ 1 GeV the allowed Sommerfeld enhancement is ~ 100. These values 
were calculated assuming a velocity dispersion of 200 km/s, and thus can be compared to our 
boost factors for the spherical halo. However, it is unclear what the results would be for the 



10 Multiple species can both decrease the needed annihilation rates for the observed PAMELA positron 
fraction by a factor of ~ 2 and evade the CMB limits. See JS^]] for a more extended discussion on indirect 
signals from multiple species of annihilating DM. 
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dark disk, which has a velocity dispersion that is ~ 5 times smaller. We note that in (55|] they 
assume the Sommerfeld enhancement reaches the saturation bound at a much lower velocity 
dispersion than that of the dark disk or the spherical halo. If the saturation bound occurs 
at velocities of ~ 200 km/s, or even ~ 30 km/s, the effect of the Sommerfeld enhancement 
on the relic density is much smaller, since much higher annihilation cross-sections at earlier 
epochs are not allowed. Thus we still consider as a reference the constraints of WMAP5 CMB 
power spectra from |99[| . We expect that data from Planck |TJ will make these constraints 
significantly tighter. 



3.4 Photons 

The predictions for the photon signals from dark matter annihilation are highly dependent on 
the distribution of the photon sources, i.e. the dark matter, in the Galaxy. Because photons 
travel to our detectors from their place of origin without diffusion and without loss of energy, 
their fluxes sample their source distribution very well. Thus, not only is the relative density 
of dark matter in the disk to that in the spherical halo (po DD / POsh ) important for calculating 
photon fluxes, but so is the relative annihilation cross-section. Whether or not there is a 
Sommerfeld enhancement and whether or not it is in saturation is quite relevant for photons. 
For a discussion of the range of reasonable values of the ratio of densities, see Section |2[ There 
are two limiting cases: Pdd /psh ~ 1 and Pdd /psh < 1- 



Microwave haze 

Cosmic ray electrons and positrons with GeV-scale energies will synchrotron radiate 
through their interactions with the Galactic magnetic field, resulting in microwave frequency 
radiation. The dark matter models that fit the PAMELA and Fermi electronic data predict 
a population of e ± with the right spectrum and spatial distribution |5^, |6^, p9[ , assuming 
a spherical DM halo) to account for the WMAP microwave "haze" found by [^] and later 
confirmed by ^8fl . We have calculated the effect of the addition of a dark disk on the syn- 
chrotron radiation in the haze region, i.e. the average value over longitudes —10° to +10° for 
latitudes between 6° and 15° from the Galactic Plane. 

Here we assume locally equal, Pdd / PSH ~ lj and spatially uniform annihilation cross- 
sections for the DM particles in the dark disk halo and the spherical halo. This case trivially 



applies for annihilating DM models like Kaluza-Klein DM as presented in [96, |28| and studied 



by fl8| , 65, 67] (among others), where the Sommerfeld enhancement cannot be used to explain 
large boost factors. As we have already discussed, the Sommerfeld enhancement may be at or 
very close to saturation, so the assumption of equal annihilation cross-sections may be valid 
for theories that include the Sommerfeld enhancement as well. However, the annihilation 
cross-section may drop near the Galactic Center, since recent simulations in which baryons 



are included [^J, |93|, 87, suggest that the velocity dispersion grows toward the GC. [87] find 
that the dependence of the velocity dispersion a v on the distance from the Galactic Center r 

-1/4 



?oes as <7„ ~ r 
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For the dark disk and Einasto profiles of Eqs. 2.1 and 2.2, under the assumption of equal 
local densities, the flux of synchrotron radiation for the spherical profile is ~ 20 times larger 
than that for the dark disk within the region of the microwave haze. This value is nearly 
constant over the entire region. When averaging over the inner 2 kpc, the Einasto profile 
has a number density squared value n 2 that is ~ 60 times larger than that of the dark disk 
60, and about ~ 5 times larger than the mixed term (see Eq. 2.4), 
5. If we define p' SH to represent the DM density in the case where 



profile, n 2 SH /n 2 DD 



Psh/ 2 Psh ■ PDD 

there is only a spherical halo, then averaging over the inner 2 kpc, assuming p' SH = p$H and 
POsh ~ Po DD 7 we nn d (psh + Pdd) 2 I p'sh ~ 5/4. Since we measure the total local DM density, 
we need to normalize the total DM density locally to the same value, po SH + po DD = Pq sh - 
For p s H ~ Pdd we get, f/ Qgg = 2p 0sH -» p' SH = 2p S u- Thus after averaging over the inner 
2 kpc (psh + Pdd) 2 / p'sh ~ V^- ^ we a ^ so normalize the BF to fit the local fluxes and 
include the effects of propagation of the highest energy e ± , then in non-Sommerfeld cases 
and saturated Sommerfeld cases we still need approximately the same BF to explain the 
WMAP haze, even including a dark disk. However, if the Sommerfeld enhancement for the 
disk is significantly greater than for the spherical halo (and still within constraints), then the 
synchrotron radiation from the e ± of DM origin can be decreased by a factor of ~ 2. 11 

Hi: 



Common variations of dark disk profile of Eq. 2.1 include 



p(R,z) 



cxp 



R 



R 



1/2 



exp 



Zl/2 



and p(R, z) ~ exp 



R 



R 



1/2 



sech 



2l/2 



(3.3) 



Taking z 1 / 2 = 1.5 kpc and R\ii = 11.7 kpc, we calculate synchrotron fluxes similar to those 
for our reference dark disk. Thus, the above arguments apply for these variations as well. 
p£[ | has suggested a dark disk with zy 2 = 2.8 kpc and R\/2 = 12.6 kpc, using the HI gas 
distribution as an alternative tracer of the gravitational field. Using these values of the scale 
lengths in our profile of Eq. 2A results in an increase in the synchrotron flux from the dark 
disk of 50% over the flux calculated using the values of Z\i 2 = 1-5 kpc and R\i 2 = H-7 kpc. 
Thus, the specifics of the profile of the dark disk have a small influence on our conclusions. 
We note that since the Fermi (7-ray) haze [49] is the 7-ray counterpart of the WMAP haze, 
the previous discussion applies to the Fermi haze as well. 



Diffuse 'j-rays 

In addition to producing synchrotron radiation of microwave wavelengths through inter- 
actions with the Galactic magnetic field, the population of e^ from dark matter annihilation 
can produce gamma rays through inverse Compton scattering (ICS) off of the optical, infrared, 
and CMB photons in the interstellar medium. Another source of 7's from DM annihilation 
is prompt emission, gammas produced in the decay chains of the annihilation products. 

n If for the dark disk the annihilation cross-section is close to a resonance, then the effect on the synchrotron 
radiation from the e of DM origin is even more prominent. 
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Figure 10: 7-ray flux from ICS of DM e ± off the ISRF as a function of galactic latitude for the XDM 
e ± channel with m x — 1.2 TeV. Left: E 7 — 10 GeV, saturated Sommerfeld enhancement. Right: 
Sommerfeld enhancement with annihilation in the dark disk only and in the spherical DM halo only 
for E 1 = 10 and 100 GeV. The isotropic diffuse 7-ray flux measured by Fermi [|| at 10 GeV and 100 
GeV(extrapolated) is shown within 2c 



In Fig. |lO| a we show the flux of 10 GeV 12 ICS 7-rays for the XDM channel, averaged 
over longitudes \l\ < 5°, as a function of latitude b. Here we assume equal local DM densities 
in the spherical and disk halos, and normalize the local e + + e~ flux to fit the Fermi data. We 
also indicate in the plot the isotropic diffuse 7-ray flux at 10 GeV as measured by Fermi ||] 
with 2c7 errorbars. 13 As recently suggested by [ 43 1 , about ~ 20% of the isotropic diffuse 7-ray 
flux for E > 0.1 GeV is attributed to unresolved extragalactic sources, mainly blazars, with a 
power-law of E~ 2 - 2 . (See also ||.) Thus, it is possible that some of the isotropic diffuse 7-ray 



flux is coming from dark matter annihilation or decay. For the scenario shown in Fig. 10a, 
the galactic diffuse 7-ray flux for \b\ > 70° is very similar in magnitude to the isotropic diffuse 
flux. It is thus evident that constraints can be put on the annihilation rate of DM models 



from diffuse 7-rays at high latitude. On a related note, |104| and j| have suggested that dark 
matter substructure at z > can contribute an important fraction to the isotropic 7-ray flux. 

have shown that calculations of the extragalactic 7-ray flux from cosmological DM 



can also put constraints on the annihilation cross-section of various DM models. 

Including a dark disk increases the galactic diffuse 7-ray flux at high latitudes (|/| > 60°) 
by 1/2 relative to the flux at low latitudes (near the GC), in cases where the main production 
mechanism of 10 GeV 7-rays from dark matter annihilation is ICS of e ± off the Interstellar 



Radiation Field (ISRF). This is readily seen in Fig. 10a by comparing the flux from the 



Spherical Halo to the total flux; the flux from the SH falls off more quickly with latitude. 



12 E*y ~ 100 GeV is also relevant for probing DM, but Fermi statistics are very low and CR contamination 
is more significant around 100 GeV. 



13 The flux is based on the data shown in 0] and 
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Therefore, a significant disk component of dark matter could result in stronger constraints 
coming from the isotropic diffuse 7-ray flux than from the diffuse flux at the Galactic Center. 
For DM models producing a large fraction their total gamma rays at EL = 10 GeV coming 
from prompt emission, the relative increase of high latitude flux to low latitude flux is smaller. 
We therefore conclude that a dark disk can't be a significant issue in the search for 7-ray 
anisotropies due to DM clumps [j47|, [75|, |97|, |4], |74|, [7(], |98|, at high latitudes. 

We now consider a situation in which the Sommerfeld enhancement is not in saturation, 
so that the annihilation cross-sections in the disk and spherical halo may be quite different. 
We assume Po DD / Po SH ~ 1- In Fig. lCb, we show the 7-ray flux from ICS of DM e^ off the 
ISRF as a function of galactic latitude for the two extreme cases, one in which the annihilation 
in the spherical halo is dominant, i.e. we ignore any dark disk contribution, and one in which 
the annihilation in the dark disk is dominant, i.e. we ignore any spherical halo contribution. 
We show results for the XDM annihilation channel with M x = 1.2 TeV, and also include 
the Fermi isotropic flux at Ej = 10 GeV and 100 GeV. We normalize the annihilation rates in 
both cases to give agreement with the local e + e~ CR spectra. The diffuse 7-ray flux at high 
latitudes is larger by a factor of 2(4) at E~ = 10(100) GeV in the disk-dominated scenario 
than in the spherical halo-dominated scenario. Thus, even if the Sommerfeld enhancement is 
currently not saturated, the existence of a dark disk with a velocity dispersion of ~ 30 km/s 
does not significantly increase the high latitude diffuse 7-ray flux. Clearly, Fig. [lO| b indicates 
that the annihilation rate in the dark disk is highly constrained by the isotropic diffuse flux 
at £ 7 = 10 GeV. 

Since the velocity dispersion increases with decreasing galactocentric distance, annihila- 
tion cross-sections toward the center of the Galaxy will be smaller in models of DM with a 
Sommerfeld enhancement, provided the enhancement is not saturated. So, the constraints on 
the BF needed to fit the local CR fluxes coming from the 7-rays measured by HESS in the 
Galactic Ridge (GR) region and inner 0.1° of the Galactic Center [12, 11, 13, 16] are lifted, as 
was shown in ]33| . While the existence of a dark disk boosts the DM signals of local origin, 
the signals, while at the same time not dramatically increasing the 7 signals from the GC, 
the real gain to be had in Sommerfeld-enhanced models comes from the effects of the increase 
in velocity dispersion toward the GC. For DM models without Sommerfeld enhancement, the 
constraints placed on the local annihilation rates by the gamma ray fluxes in the GR and GC 
regions [20, [22|, |36j, |45l] are similar with and without a dark disk component. 



4. Conclusions 



A self-consistent explanation of the PAMELA and Fermi signals as coming from dark matter 
annihilation requires a local annihilation rate in the spherical halo that is 2 to 3 orders of 
magnitude larger than that expected for a thermal WIMP. The large boost factor could arise 
from substructure in the halo (a local inhomogeneity) , from the Sommerfeld enhancement, 
or from a Breit-Wigner resonance, or it could be an indication of a non-thermal WIMP. We 
consider effects of a dark disk component with a local density of 0.4 GeVcm -3 . The lower 
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velocity dispersion of dark matter particles in the dark disk can naturally lead to higher 
annihilation rates through the Sommerfeld enhancement, which scales as l/v re i. We show 
that if the dominant dark matter contribution to the local high energy e^ flux comes from a 
dark disk, rather than the spherical halo, we can get good fits to the Fermi and PAMELA with 
a boost factor that is almost an order of magnitude smaller than what is needed for the fits 
to the data when only the spherical halo profile is taken into account. Thus, the Sommerfeld 
enhancement for the spherical halo with a velocity dispersion of v re \ ~ 200 km/s need not take 
the extremely large values previously thought. We stress that many of our conclusions would 
hold to varying degrees for all flattened halo profiles, even those with kinematics similar to 
the spherical halo, since all of the cosmic ray dark matter signals are highly dependent on 
the spatial distribution of the dark matter. 

We show that the locally observed spectra of e + e~ are harder for a dark disk profile than 
for a spherical profile. This can improve the fits to the Fermi data for those annihilation 
modes with softer e + e~ spectra. Some tension arises as a result of adding a disk component 
to the dark matter density, because the discrepancy between the boost factors needed to fit 
the PAMELA positron fraction and to fit the Fermi data increases with the addition of a disk 
component. The spherical halo gives differences of factors of ~ 2 — 3, while the dark disk gives 
differences of factors of up to ~ 9 (see Figs. One explanation is that the spherical halo 

contribution to the e + e~ flux is dominant at the lower energies measured by PAMELA, while 
for energies greater than 300 GeV, the Fermi range, the dark disk contribution to the e + e~ 
flux becomes comparable. That, however, depends on the assumptions about the contribution 
to the local DM density from the dark disk and the relative annihilation cross-sections in the 
dark disk and the spherical halo. Furthermore, the possibility of multi-component DM or 
nearby pulsars contributing significantly to the local spectra decreases the significance of a 
difference in the boost factors necessary to fit PAMELA vs Fermi. 

We find that for hadronic annihilation channels, if the annihilation in the disk is the 
dominant contribution to the local cosmic ray fluxes of dark matter origin (and if the disk 
is thicker than the diffusion zone), the data allows for an annihilation rate that is larger by 
a factor of ~ 6 than the rate for a spherical profile only. Even so, the PAMELA antiproton 
data constrains the annihilation products to be mainly leptonic for masses of m x ~ 1.0 TeV. 

We investigate possible constraints on the dark disk contribution to the e^ fluxes using 
the Fermi e + + e~ data and CMB data. The Fermi data is consistent with a scenario in which 
the dark disk is the dominant contributor to the local e + e~ flux for the XDM annihilation 
channels we studied. With the data currently available, it is impossible to determine whether 
the local e + e~ flux is produced dominantly in one halo component or approximately equally 
in both the disk and the spherical halo. In fact, due to the many possible dark matter masses 
and annihilation modes in the production of the local e + e~, it will be very hard to tell the 
difference even with better data. The constraints on the DM annihilation rate coming from 
the WMAP5 measurements of the temperature and polarization power spectra of the CMB 
do not exclude any of the XDM annihilation channels studied in a disk-dominated scenario. 
However, because the dark disk produces fewer electrons and positrons at PAMELA energies 
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than the spherical halo, there is more tension for the disk than for the spherical halo between 
the BF needed to fit the PAMELA data and the CMB constraint. 

We find that the contribution from the dark disk can be up to ~ 1/2 of the local 
flux of dark matter origin for E e ± > 500 GeV. Additionally, we find that including a dark 
disk has a minimal effect on the flux of synchrotron radiation from dark matter in the 
microwave haze region. This result is valid even for disks with relatively large scale heights 
(~ 3 kpc) perpendicular to the plane. Moreover, including a dark disk of scale height ~ 1 kpc 
does not affect the diffuse 7-ray flux at high latitudes by more than 0(1). If the Sommerfeld 
enhancement is not saturated, so that the annihilation rate in the disk is much greater than 
the annihilation rate in the spherical halo, the contribution to the gamma ray flux from the 
disk at 10 GeV exceeds the value of the isotropic 7-ray flux as measured by Fermi by a factor 
of ~ 2 at the highest latitudes where the isotropic component is the dominant component. 
Thus, the annihilation rate in the dark disk is highly constrained by the isotropic diffuse flux 
at E 1 = 10 GeV. 
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The probability distribution function of injected electrons created by the annihilation 
channel \X ~~ ^ ~~ * 4</> — >■ (4/i 1 * 1 —¥) 4e ± can be calculated analytically or numerically. In 
the frame of (ft, it is the convolution of the spectra of electrons and positrons in that frame with 
the spectrum of (ft as seen from the frame in which annihilation takes place. For simplicity in 
our analytic calculations, we consider <E> and (ft to be scalars so that the decay products have an 
isotropic distribution in the rest frame of the parent particle. The difference in the spectra 
resulting from a vector mediator is a spectral hardening of O(0. 1) at the highest ~ 10% of the 
energies. Since the deviations appear only at the highest energies where experiments have 
their largest errorbars, the spectra for the scalar mediator case can be considered a good 
approximation to the spectra for the vector mediator case. In both cascades described, the 
multiplicity of electrons (positrons) produced per annihilation is 4. Our results agree with 
those of Q. 

For the "direct" case, (ft — > e^, the distribution in the approximation m x ^> m$ ^> ^ 
m e is given by: 



Appendix A: 2-Step e + e spectra 




(1) 
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For the case of <p—> and consequently muon decay, the spectrum's given by 
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